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Quantum mechanics and classical mechanics are two very different theories, but the correspon- 
dence principle states that quantum particles behave classically in the limit of high quantum number. 
In recent years much research has been done on extending both quantum mechanics and classical 
mechanics into the complex domain. This letter shows that these complex extensions continue to 
exhibit a correspondence, and that this correspondence becomes more pronounced in the complex 
domain. The association between complex quantum mechanics and complex classical mechanics is 
subtle and demonstrating this relationship requires the use of asymptotics beyond all orders. 
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The correspondence principle states that at high en- 
ergy quantum mechanics resembles classical mechanics. 
In Fig. [I] we illustrate this resemblance for the harmonic 
oscillator Hamiltonian H — p 2 + x 2 . We compare the 
quantum probability density p qU antum 0) = \ipi 6 (x)\ 2 for 
the 16th energy level E\q = 33 with the classical prob- 
ability density p c iassicai(^) = wnere s is the speed 
of the particle. In the classically allowed region, which 
is bounded by the classical turning points at ±\/33, 
PquantumO) oscillates about PciassicaiO)- The particle 
speed vanishes at the turning points, so ^classical is singu- 
lar. (This singularity is integrable and thus the classical 
probability is normalizable.) In the classically forbidden 
region p qU antum(^) decays exponentially; in conventional 
classical mechanics the particle does not enter this region 
and its probability density is thought to vanish. 

This paper generalizes quantum and classical proba- 
bility into the complex domain. In complex classical me- 
chanics we solve Hamilton's equations x = p = — ^ 
for complex initial conditions and not just for initial con- 
ditions in the classically allowed region [TJ [2j |3j HJ El El [7] . 
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FIG. 1: Correspondence principle for the harmonic oscillator. 
The probability density p(x) — |-0i6(x)| 2 for a quantum par- 
ticle (solid curve) and the corresponding probability density 
for a classical particle of the same energy (dotted curve) are 
shown. In the parabolic potential well p qU antum(^) is wave- 
like and oscillates closely about p c iassicai(^). At the classical 
turning points, /^classical (^) diverges. In the classically forbid- 
den region the quantum particle density decays exponentially 
while the classical density is ordinarily assumed to vanish. 




FIG. 2: Classical probability density at z — x + iy in the 
complex plane for a particle of E — 2 in a harmonic potential. 
The probability density resembles a pup tent with infinitely 
high tent poles located at the turning points. The classical 
trajectories are nested ellipses, all of which have the same pe- 
riod T = TV. These complex orbits are superposed on the tent 
canopy. The degenerate ellipse, whose foci are the turning 
points at x = ±y/~E, is the conventional real solution. Clas- 
sical particles repeatedly cross the real axis in the classically 
forbidden regions |x| > y/~E. 



For the harmonic oscillator the classical orbits are nested 
ellipses with foci at the turning points. Using p qU antum in 
Fig. [T] as a guide, we require the probability of a classical 
particle being on more distant ellipses to decay exponen- 
tially, and we plot in Fig. [2] the relative probability den- 
sity of finding the classical particle at the point z = x+iy. 

Figure [2] shows that, contrary to the traditional view, 
the classical particle density in the classically forbidden 
region is actually nonzero, and that ^classical resembles 
Pquantum in this region [8 . This figure emphasizes that 
the classical probability density extends beyond the real 
axis and into the complex plane, where p c iassicai(^? y) falls 
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off as the reciprocal of the distance from the origin. 

To establish the complex correspondence principle we 
must continue the quantum probability density into the 
complex plane to match the complex classical probabil- 
ity density in Fig. [2j To do so we must extend quan- 
tum mechanics into the complex domain [9 . A Her- 
mitian Hamiltonian (H = H^) has real energy levels 
and unitary time evolution. However, the class of phys- 
ically allowable Hamiltonians may be extended to in- 
clude non-Hermitian Hamiltonians that possess an un- 
broken VT (combined parity and time-reversal) symme- 
try [101 HQ [I! EH HI US] - These complex Hamiltonians 
also have real spectra and generate unitary time evolu- 
tion, and such Hamiltonians have recently been observed 
in the laboratory [T6l fT7| H5] . 

For T^T-symmetric Hamiltonians the potential satisfies 
V*(— z) = V(z). This allows us to derive a conserva- 
tion law in the complex plane from the time-dependent 
Schrodinger equation iip t (z,t) = —ijj zz (z,t)+V(z)ilj(z,t). 
The continuity equation is pt(z,t) + j z {z,t) = 0, where 
the local probability density and current are 

p(z,t) = i/>*(-z,t)rl>(z,t), (1) 
j(z,t) = (2) 

In this paper we limit our discussion to wave functions 
ip that are eigenstates of H; for such states the density 
p(z) is time independent and the current j(z, t) vanishes. 
The density p(z) is not the absolute square of ip(z, t) and 
it is analytic in the complex- z plane. 

A locally conserved probability density must be real 
and positive and its spatial integral must be normalized 
to unity. While p(z) satisfies a continuity equation, it 
cannot be a probability density because it is complex- 
valued. We propose a novel approach. We construct a 
complex contour C satisfying three conditions: 



Condition I : 
Condition II : 



Im [p(z) dz] = 0, 
Re \p(z) dz] > 0, 



Condition 



III : / Re[p(z)dz] = 1. 
Jc 



(3) 
(4) 

(5) 



For brevity, we discuss the harmonic-oscillator ground 
state tpo(z) = e~ z I 2 for which p(z) = e~ x +y ~ 2lx y . 
Since dz = dx + idy, Condition I gives 



dy/dx = $m(2xy)/ cos(2xy), 



(6) 



which is a differential equation for the contour C. On 
the contour defined by this differential equation the local 
contribution p(z) dz to the probability is real. 

Next, we turn to Condition III. Inserting ([6| into (J5|, 
we get two forms for the probability integral: 



r e [y(x)] 2 -x 2 r 
Jc cos[2y(x)x] Jc 



e y 2 -[<y)¥ 
c dV sm[2yx{y)}' 



These integrals converge if the contour C terminates in 
the two good Stokes' wedges of opening angle tt/2 cen- 
tered about the real axis, but they diverge if C terminates 
in the bad Stokes' wedges of opening angle ir/2 centered 
about the imaginary axis. 

To determine whether the contour C terminates in 
good Stokes' wedges we must solve (J6|. However, this 
simple-looking differential equation does not possess a 
closed-form solution, and we must rely on asymptotic 
techniques. For large x the contour C approaches the 
center of the good Stokes' wedge and for large C ap- 
proaches the center of the bad Stokes' wedge: 

y(x) ~ nir/(2x) (x -> +oo, n e Z), (8) 
x(y) - (ra + 1/2) 7r/(2?/) (y — > +oo, m G Z). (9) 

This asymptotic analysis shows that the integration 
paths occur in quantized bunches labeled by the integers 
m or n. Higher-order asymptotic analysis reveals that 
these bunches are stable in the bad Stokes' wedges and 
unstable in the good Stokes' wedges. This means that 
as | y | — > oo, a path in the bad Stokes' wedge is drawn 
towards the quantized curves in (J9|, but that as \x\ — > oo 
paths veer away from the quantized curves in (J8|, which 
we call separatrices. Thus, only the isolated separatrices 
ever reach oo in the good Stokes' wedges; all other curves 
turn around and are drawn into the bad Stokes' wedges. 
This remarkable behavior is illustrated in Fig. [3] 

Figure [3] shows that the only continuous path connect- 
ing the two good Stokes' wedges follows the real axis. All 
other paths terminate in one good and one bad Stokes' 
wedge, and thus the probability integral ^ along such 
paths diverges. This instability of contours in the good 
Stokes' wedge in Fig. [3] is a serious and generic problem 
that must be overcome if we wish to extend the quantum 
probability density into the complex plane. 

To overcome this apparently insurmountable problem, 
we show that if a contour satisfying ([6| in the complex-z 
plane enters a bad Stokes' wedge along one path and then 
returns along a second path in the same quantized bundle, 
then the integral along the combined paths is actually 
convergent! This result is surprising because the integral 
along each path separately is exponentially divergent. To 
prove convergence, we examine the integral 
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dy e y 



' e -[<y)] 2 e -[v(y)} 2 \ 

sm[2yu(y)} s'm[2yv(y)] J 



(10) 



(7) 



where u(y) and v(y) are two solutions to ([6| in the same 
bunch (that is, having the same value of m). The leading 
asymptotic behaviors of these solutions are given in ([9|, 
and in fact the entire Poincare asymptotic expansions of 
u{y) and v(y) are identical to all orders in powers ofl/y. 
However, the difference D(y) = u{y) — v(y) is nonzero 
and is exponentially small 

D(y) - Ce~ y2 [1 + (n + 1/2) 2 tt V 2 /4 + . . .] (11) 
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FIG. 3: Numerical solutions (solid curves) to the differen- 
tial equation Q in the complex- z — x + iy plane. Solutions 
have vanishing slope on lightly dotted hyperbolas and infi- 
nite slope on heavily dotted hyperbolas. In the good Stokes' 
wedges (dark shading) p(z) decays exponentially and in the 
bad Stokes' wedges (unshaded) p(z) grows exponentially. The 
only curves that reach ±00 in the good wedges are separ at rices 
(heavy dashed curves labeled n = 0, . . . , 4). A typical solu- 
tion curve in a good Stokes' wedge is unstable; as x — > 00 the 
curve turns away from the separatrix, leaves the wedge, and is 
drawn into a bad Stokes' wedge. In the bad wedge curves are 
stable and continue on to ±ioc in quantized bunches labeled 
by m. The only continuous unbroken curve connecting the 
two good wedges is the special separatrix on the real axis. 



as I y I — > 00, where C is an arbitrary constant. This 
result reflects the hyperasymptotic (asymptotics beyond 
all orders) content of |6| [19] . 

Using ( pTj ), we approximate the integral (10): 



7~C(-l) n (n + 



1/2) / 

JY 



dyy- 3 (y-oo), (12) 



which is convergent. Thus, while there is no continuous 
path running between the two good Stokes' wedges, there 
do exist paths connecting these wedges that repeatedly 
enter and re-emerge from bad Stokes' wedges. On these 
contours the probability integral / in ^ exists. Several 
such paths are shown in Fig. [4] 

Figure[5]is a generalization of Fig.[4]for the first excited 
state of the quantum harmonic oscillator. Note that con- 
tours that connect the two good Stokes' wedges may or 
may not pass through the node at the origin. 

We do not present the argument in this paper, but 
along the contours shown in Figs. [4] and [5] Condition II in 
Q is satisfied. That is, along the entire complex contour 
obeying the differential equation |6| the local contribu- 
tion to the total probability integral is positive. Thus, 
the requirements in the three conditions ([3]-[5| are met, 
and we have successfully extended the probabilistic inter- 
pretation of quantum mechanics into the complex plane. 




FIG. 4: Complex contours connecting good Stokes' wedges. 
A contour (solid curve) begins at x = —00 in the left good 
Stokes' wedge (shaded), leaves the wedge along a separatrix, 
and runs up to ioc in the bad Stokes' wedge (unshaded). The 
contour continues downward along a path in the same bunch, 
crosses the imaginary axis at y = 0.003 324 973 872 707912, 
and heads upwards into the same bad Stokes' wedge. Finally, 
the contour re-emerges from the bad Stokes' wedge and con- 
tinues towards x = +00 along a separatrix in the right good 
Stokes' wedge. Contours have zero slope on lightly dotted 
lines and infinite slope on heavily dotted lines. A second con- 
tour (dashed curve) connecting the two good Stokes' wedges 
visits the bad Stokes' wedge four times instead of twice. 




FIG. 5: Complex contours for the first excited state of the 
quantum harmonic oscillator. Four contours (solid lines) that 
begin at x = —00 in the left good Stokes' wedge (shaded) 
are shown. These contours leave this wedge along separatri- 
ces and run off to zbzoo in the upper and lower bad Stokes' 
wedges (unshaded). After visiting a bad Stokes' wedge several 
times, the contours may pass through the node at the origin 
at 60° angles to the horizontal. At this node the probability 
density vanishes. Then the contours repeat the process in the 
right-half plane and eventually enter the right good Stokes' 
wedge along separatrices. It is also possible to have a complex 
contour that avoids the node at the origin and still connects 
the good Stokes' wedges. Solution curves are horizontal on 
lightly dotted lines and vertical on heavily dotted lines. 
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FIG. 6: Complex contours (heavy dashed lines) for the ground 
state (J = 1) of the quasi-exactly-solvable VT anharmonic 
oscillator. A separatrix (a) goes directly from the left good 
Stokes' wedge to the right good Stokes' wedge, crossing the 
imaginary axis at y = -0.176 651795 619 462 368. Paths (b) 
that cross the imaginary axis at higher points than the (a) 
path cannot reach oo in the good wedge. Because they are 
unstable, such paths turn around, enter the upper bad Stokes' 
wedges, and can never re-emerge. A separatrix (c) leaves the 
left good wedge. This path enters the lower bad wedge to 
the left of the imaginary axis, re-emerges along paths (d) or 
(h), and reenters the bad wedge to the right of the imaginary 
axis. It then continues into the right good wedge along the 
separatrix (e). Another separatrix (f) leaves the left good 
Stokes' wedge, enters the lower bad Stokes' wedge, leaves and 
returns along (g), leaves, reenters again along (d) or (h), leaves 
and reenters along (j), and finally enters the right good wedge 
along the separatrix (f). Solution paths are horizontal on 
lightly dotted lines and vertical on heavily dotted lines. 

The analysis in this paper is general and extends to 
other VT quantum theories, such as the quasi-exact ly 
solvable quartic anharmonic oscillator, whose Hamilto- 
nian is H = p 2 - x A + 2ix 3 + x 2 + 2iJx [20]. However, 
for quantum theories other than the harmonic oscillator, 
there are two kinds of bad Stokes' wedges: (1) wedges for 
which there exist probability contours C that enter and 
re-emerge from the wedge and for which the probability 
integral along C converges and (2) wedges for which no 
such contour exists. For the Hamiltonian above there are 
three good Stokes' wedges and three bad Stokes' wedges, 
one of type (1) and two of type (2) (see Fig. |6|. 

The quantum probability densities associated with the 
curves in Figs. |1J |5j and [6] resemble pup tents with rip- 
ples in their canopies, in contrast to the smooth surface 
in the classical pup tent in Fig. [2] This is the complex 
analog of Fig. [T] and illustrates the complex correspon- 



dence principle. 
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